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We argue that the chemical freeze-out in heavy ion collisions at high baryon density is not as- 
sociated to a phase transition or rapid crossover. We employ the linear nucleon-meson model with 
parameters fixed by the zero-temperature properties of nuclear matter close to the liquid-gas quan- 
tum phase transition. For the parameter region of interest this yields a reliable picture of the 
thermodynamic and chiral properties at non-zero temperature. The chemical freeze-out observed in 
low-energy experiments occurs when baryon densities fall below a critical value of about 15 percent 
of nuclear density. This region in the phase diagram is far away from any phase transition or rapid 
crossover. 



Relativistic heavy ion collisions are a promising way to 
investigate the properties of fundamental quantum field 
theory - more specific QCD - at nonzero temperature 
and density. However, in contrast to most condensed 
matter systems the produced matter has only a short live 
time during which it expands and cools quickly. It is a 
big challenge to determine the properties of the produced 
matter. Basically all information has to be reconstructed 
from the final state, i. e. from the momenta and chemical 
composition of the detected particles. 

Important observables are the yields of different parti- 
cle species. To surprisingly good approximation they can 
be described by the so called statistical model which as- 
sumes a thermal distribution of a non-interacting hadron 
resonance gas [IH3^ • The most prominent observables ex- 
tracted from the thermal fits are the temperature and 
baryon chemical potential associated with the chemical 
freeze-out. 

It has been advocated that the chemical freeze-out 
temperature coincides with the temperature of the QCD 
phase transition for low baryon density |3] . The basic ar- 
gument states that the rates of particle number changes 
in the hadronic phase are too small to maintain chem- 
ical equilibrium. This holds whenever scattering pro- 
cesses involving only a few particles dominate. Chemical 
freeze-out therefore occurs for particle number densities 
that are just high enough that multi-particle scattering or 
collective effects dominate. For low baryon number that 
is possible only for temperatures very close to a phase 
transition or crossover. The difference between the criti- 
cal temperature and the freeze-out temperature has been 
estimated to be less than 5 MeV. The measured freeze- 
out temperature at RHIC at heavy ion collisions with 
= 200 GeV of 164 ± 6 MeV [3] agrees actually 
rather well with the critical temperature of the crossover 
in lattice simulations of Tc « (157 ± 10) MeV [5] and 
Tc « (154 ±9) MeV [g. 

The agreement between the critical temperature of a 
phase transition or crossover and the observed chemical 
freeze-out temperature is expected to hold for low values 
of the chemical potential fi in the QCD phase diagram. 
Naturally, the question arises whether a similar argument 
can be extended to higher values of ^ or larger baryon 
density. Interesting ideas for this issue have been pro- 



posed recently [7^. We ask: does the curve of measured 
freeze-out temperatures reflect a phase transition line or 
rapid crossover in the whole //-T-plane? 

In this letter we argue that this is not the case. The 
observed freeze-out temperatures T^h for the largest val- 
ues of ^ lie actually in a region where a simple modeling 
by baryons and mesons becomes possible. While a rapid 
change of the particle density with temperature continues 
to play a crucial role for the determination of T^h, this 
actually happens in a region that is substantially away 
from any transition or crossover. We illustrate the situa- 
tion in Fig.[l]where we indicate the observed points in the 
/i-T-plane and the region of validity of a simple baryon- 
meson model. We also demonstrate in Fig. [2] the change 
of particle density with temperature for fixed value of /x. 
The dot in this figure indicates the measured value of Tch 
for fi — 760 MeV. A similar figure [3] for the chiral order 
parameter ctq as a function of T shows that no particu- 
lar distinct feature such as a (chiral) phase transition or 
crossover is visible in this range. The deviations of <tq 
from the vacuum value are small in the whole range of 
the black solid curve in Fig. [l] 

For a discussion of the phase diagram and the thermo- 
dynamic and chiral properties within the region indicated 
in Fig. [TJthe linear nucleon-meson model is a reasonable 
approximation. It uses as degrees of freedom the pro- 
ton and neutron, a neutral vector meson o;^ as well as 
the pions and the collective cr-meson. (We ignore isospin 
violation and electromagnetism for simplicity.) Chiral 
symmetry is implemented explicitly. Integrating out the 
cr-degree of freedom yields an effective non- linear ct- model 
coupled to nucleons and the vector meson. For this lat- 
ter model the chiral perturbation theory has been used 
extensively [S]. On the other hand, within a quadratic ap- 
proximation to the effective potential of the field cr, and 
if only the dominant nucleon fluctuations are included 
for the computation of the chiral order parameter and 
the baryon density, one recovers the gap equations of the 
Walecka model [S; . Parameters of the effective potential 
at zero temperature and the chemical potential /i^ cor- 
responding to the gas-liquid phase transition in nuclear 
matter can be determined from observation. In this pa- 
rameter region the model can be mapped directly to the 
nuclear droplet model. 
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FIG. 1: Curve of constant baryon number riBaryons = 
0.15 Jinucicar lu the Mesou-Baryon model (solid black line). 
The points with error-bars mark the chemical freeze-out as 
obtained from the fits to experimentally measured particle 
yields [5]. The red line marks the first order phase transition 
to nuclear matter. The dashed and dashed-dotted lines indi- 
cate an estimate for the range of applicability of our model. 
More specific, in the region to the right of the dashed line the 
relative contribution of pions to the pressure is smaller than 
20%. In the region to the left of the dashed-dotted line the 
baryon density nearyons is smaller than 1.5 times the nuclear 
saturation density nnuciear = 0.153/fm^. In this region no 
signs of a phase transition are visible. 
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FIG. 2: Number density of baryons as a function of the 
temperature for /i — 750 MeV (solid line). Note that the 
number of anti-baryons is negligible within the plot resolu- 
tion. We also show the number of pions (dashed line). The 
dot marks the experimental result for the chemical freeze- 
out temperature Tch = 56^.^'^ MeV corresponding to Hch = 
760 ±22.8 MeV. 



The computational task concerns then mainly the dif- 
ference of the effective meson potential U{a;T,fi) — 
U{(j;0, fic)- This can be done by various methods - for 
example one could employ functional renormalization by 
adding nucleon degrees of freedom to the setting of ref. 
[lOj . For our limited purpose a very simple approach 
will do. The potential difference is directly related to 
difference of pressure for the parameters {a; T, /i) and 
{a;0,fj,c). This can be approximated by a free gas of 
nucleons with cr-dependent mass. We can consider a as 



FIG. 3: Chiral order parameter as a function of the tempera- 
ture for jj. — 750 MeV. The dot marks the experimental result 
for the chemical freeze-out temperature Tch = 56 
corresponding to /ich = 760 ± 22.8 MeV. 
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an additional parameter in thermodynamics. Its value 
can be varied by varying the quark mass. If needed, me- 
son fluctuations can be added in a similar way. We will 
discuss the linear nucleon-meson model in the setting of 
ref. [TT]. (Our normalization of a differs by a factor 2 
from |lll.) Our new results extend the analysis to non- 
vanishing temperature. 



Linear nucleon-meson model 

We use an effective model for baryons tpa (a is an 
isospin index with ipi describing protons and '02 neu- 
trons), an isospin singlet vector meson w^, a scalar meson 
a and pseudo-scalar mesons tt^ = tts, tt* — ^(7ri±i7r2). 
It is convenient to combine the scalars and pseudo-scalars 
in the field 



L(a-i7rO) 



(1) 



V2 



The effective Lagrangian is of the form 

+ V2h[i>a (^) Cl^abi^b + (^) ('/'^)abV'6] 

-|-i0:,(-a^9^)0afc-f [/^ic(p,a) 



(2) 



Here we use the chiral invariant scalar field combination 
p — ^(/'*f,'/>ah and C/mic(p, cr) is a microscopic form of the 
effective potential 



Urnicip, o") = U{p) ~ mlUa. 



(3) 



The Lagrangian (pi is invariant under the chiral symme- 
try SU{2)v X SU(2)a X U{l)v x U{1)a where the nucleon 
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doublet transforms according to 



11 It 



It % % 



(4) 



where t — {ti,T2,t^) are the Pauli matrices in isospin 
space and the scalar field transforms according to 

(l> -)■(/>- ^av[T,(f>\- -aA{T,(l)} + il3A(t>. (5) 

The vector meson field is invariant under the above 
symmetry. It couples to the conserved baryon number 
current associated with the U{l)v symmetry such that 
only the spin-one part of plays a role while the spin- 
zero component S^w^ decouples. 

The only explicit breaking of chiral symmetry comes 
from the quark masses. This is reflected by the linear 
term in the effective potential Therefore f/mic de- 
pends explicitly on the field a in addition to the invariant 
p = ^(cr^ -I- TT^). The scalar field a will have a vacuum 
expectation value, in contrast to the pseudo-scalar field, 
TT^ = TT"*" = 7r~ = 0. Due to rotational symmetry, only 
the zero-component of the vector field can have an 
expectation value. Inspection of Eq. ^ shows that this 
can be interpreted as a shift in the effective chemical po- 
tential. 

We are interested in the quantum effective potential 
f/((T, Wo) which includes the effects of quantum and ther- 
mal fluctuations. It can be obtained from the quan- 
tum effective action - the generating functional for one- 
particle irreducible Greens functions - by specializing to 
constant a and uj with — 0,?: — 0. The minimum of 
?7((T, ujq) determines the expectation values for a and ljq, 
i.e. the chiral order parameter and the effective chemical 
potential. In general, the computation of U (cr, wq) from a 
microscopic action is a complicated task. From the real- 
ization of symmetries we know, however, that the explicit 
symmetry breaking occurs only via the unaffected linear 
term, 

U{(T,ujo) = U{p,ujn)-mlf^a , p=^a'^, (6) 

such that the task consists in a computation of the chi- 
rally invariant potential C/(p, ojo). 

We are only interested in the difference A = 
U(p,ujq;T, p.) — U{p,u!o;0, pc), with pc the value of the 
chemical potential at which the zero-temperature phase 
transition between a hadron gas and nuclear matter oc- 
curs. This simplifies our task considerably. Instead of 
a complicated computation of U{p,uiQ]Q, pc) we can use 
observation in order to pin down the relevant proper- 
ties of this quantity. In a functional renormalization ap- 
proach the difference A involves only mesons with mass 
m smaller ttT or baryons with m — pcS smaller than ttT, 
where p^f f = p + gu)Q and Qq is the expectation value of 
Wo [ini EJ • In our range of interest these are essentially 



nucleons and possibly pions, thus justifying the degrees 
of freedom incorporated in our model. 

Furthermore, for the relatively narrow range in T and 
p that we investigate here the running and the associated 
p- and T-dependence of the couplings h and g is small 
and can be neglected. In a first approach we also neglect 
the subleading pion fiuctuations. With these approxima- 
tions the solution of functional flow equations actually 
reduces to performing a Gaussian functional integral over 
the fermionic fields in a background of constant a and 
ujq. This is nothing else than relativistic mean field the- 
ory. We stress that mean field theory is generically not 
expected to give reliable results in our setting with strong 
interactions. For example, a mean field computation of 
U{p,ujq;0, pc) would fail badly. However, the more gen- 
eral view from a functional renormalization perspective 
permits to asses that the mean field result for A is reli- 
able within an appropriate parameter range. We expect 
leading corrections from the omitted pion fluctuations. 
(They can be incorporated, in principle, in some type of 
extended mean field theory, see below.) To the right of 
the left dashed line in Fig. [T] the pion contributions to 
the pressure are less than 20 %. A second type of cor- 
rection is expected due to the neglected cr-dependence of 
h, g and m^. In view of the small deviation of a from its 
vacuum value visible in Fig. |3]we expect this effect to be 
small. It increases, however, for larger density and this 
is one of the restrictions for the limitation of validity of 
our model indicated by the right dashed line in Fig. [l] 

We next discuss the mean field contribution to A which 
is directly related to the pressure of a free nucleon gas 
with field dependent masses. The corresponding contri- 
bution to the effective potential depends on the tempera- 
ture T and chemical potential p. It can be parametrized 
in terms of the pressure of a free gas of relativistic 
fermions and corresponding antiparticles 



, 1 /" cPp p 



1 



(7) 



Within our model, the effective potential for the bosonic 
fields reads now 



U{a,ujo-T,p) =J7vac(o-,a;o) 



4 Pfg{T,p + guj(),h(T), 



(8) 



where the factor 4 accounts for the degeneracy in spin 
and isospin. For the effective potential in the vacuum (i. 
e. at T = p = 0) we use the parametrization 

C/vac(a,wo) =im2(2p- /2) + h{2p ~ f^)' 

1 73 , 



+ -jl(2p-f^r + \^^i2p-f^y (9) 
1 



mlUia - U) - -mlujl. 
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For T = the pressure Pfg(0, /i, rn) vanishes identically 
for n < m. Thus A vanishes for fi < ha — gujQ and the 
expressions in Eq. ([8| and ^ coincide. 



Parameters 



At this point the open parameters of the model are 



m^, A, 73 , 74, 



h and 



(10) 



We choose the experimentally established values /tt = 
93 MeV for the pion decay constant, ni-^ = 135 MeV 
for the mass of pions and rriu: = 783 MeV for the mass 
of the vector meson. The Yukawa coupling h is fixed 
by the requirement that hf-^ equals the nucleon mass 
rUn = 939 MeV which gives h — 10. We have verified 
that somewhat smaller values of h (cf. ref. jlTj) do not 
change our results qualitatively. 

The remaining open parameters g, A, 73 and 74 are 
fixed by requiring that the model describes at vanish- 
ing temperature T = normal nuclear matter at the 
gas-liquid phase transition. More specific, the minimum 
of the effective potential Q is at cr = for vanishing 
chemical potential fj, and the baryon density hb — ^'§^U 
vanishes at this point. However, for increasing /i a sec- 
ond minimum at a smaller value of a will develop and 
a first order phase transition takes place when the two 
minima are degenerate. From the nuclear binding energy 
fbind — —16.3 MeV one can determine the critical chem- 
ical potential, ^lc = 939 MeV - 16.3 MeV = 922.7 MeV. 

From minimizing the effective potential ([s]) with re- 
spect to Wo one finds the self-consistency equation 



Wo = 



9 



(11) 



with the baryon density 



nBiO,H + gc^o,hcr) = A—Pfg{0, fi + gujQ, ha). (12) 

At fj,c the baryon density jumps fro zero to nuclear sat- 
uration density rinuci = 0.153/fm'^ this shows that also 
Wo changes discontinuously from wq = to wo.nuci at 
the first order phase transition. From the Landau mass 
niL = V^^o-nuci +Pf = + pi^o.nuci One can determine 
g and wo,nuci- We use = 0.80 m„ for our numerical 
calculation, which gives g = 9.5 and Wo.nuci = —18 MeV. 
Using the relation rinuci = ^^Pp, one can also deter- 
mine the position of the second minimum of the effec- 
tive potential with respect to a at the phase transition, 
CTnuci = 69.8 MeV. 

Of the remaining three parameters A, 73 and 74 two get 
fixed by the constraints for a first order phase transition 



J/(a'nucl,Wo,nuci;0,/lc) = [/(A , 0; 0, /ic) (13) 



and 



The third parameter (say A) can now be adapted to other 
properties of nuclear matter. For example, the choice 
A = 50, 73 = 3 and 74 = 50 corresponds to the (vacuum) 
mass of the cr-meson nic — \J A/^ -I- m\ = 670 MeV, the 
compressibility module K = 9n/{dn/dfj,) = 300 MeV 
and the surface tension of a nuclear droplet S — 
^ll^^^ y/2U{a)da = 42000 MeV^ Considered the sim- 
plicity of the model, these values are in reasonable 
agreement with the experimentally established values 
m„ « 484 ± 17 MeV [H, fsT « 240 ± 30 MeV and 
S « 42200 MeV^ 

As an independent check we determine for this choice 
of parameters also the "nuclear cr-term" which quantifies 
how the vacuum mass of the nucleon depends on the chi- 
ral symmetry breaking explicit mass term for the quarks. 
It can be determined by comparing the expectation value 
CTnuc in the nuclear matter phase with the value obtained 
in the chiral limit, = 0, for otherwise identical C/vac- 
One finds 40 MeV, in reasonable agreement with lattice 
calculations [H]. Overall, one finds a satisfactory agree- 
ment with nuclear matter at vanishing temperature and 
the nuclear droplet model. Our parameters are found in 
a similar range as the ones determined in ref. [11]. We 
have checked that other reasonable parameter choices do 
not modify our main conclusions. 



Thermodynamic properties and chemical freeze-out 

Let us now investigate the linear nucleon-meson model 
at non-zero temperature. We start from the effective 



da 



C^(o'nucl,Wo,nuci;0,/Zc) = 0. 



(14) 



potential U{a; T, /x) which follows from Eq. ( 10 ) by min- 
imizing with respect to the value of wq. The complete 
information about the phase diagram and the thermody- 
namic properties is encoded in this quantity. For exam- 
ple, the chiral condensate crp {T, /i) is determined as the 
global minimum oiU{a; T, with respect to ct. A phase 
transition of first order occurs when the potential has 
two local minima, such that at the critical temperature 
Tc the global minimum jumps discontinuously from one 
local minimum to the other. To illustrate this behavior 
we plot in Fig. |4] the effective potential U{a;T,^) as a 
function of a for vanishing temperature T = and differ- 
ent values of the chemical potential ^. One can clearly 
see the first order phase transition at the critical chemical 
potential /^^ = 922.7 MeV. 

Let us now follow the changes in the effective potential 
as the temperature is increased. In Fig. [5] we plot U[a) 
for /i = fJ.ciT) and for different temperatures T = (/ic = 
922.7 MeV), T = 5 MeV {nc = 921.3 MeV), T = 10 MeV 
{fic = 917.2 MeV), T = 15 MeV (^ic = 910.8 MeV) and 
T = 20 MeV = 902.3 MeV). One finds that the ef- 
fective potential at the minima gets more negative as the 
temperature increases. This corresponds to an increase 
in pressure p — —Umin- One also finds that the poten- 
tial barrier between the two minima becomes smaller, 
such that the surface tension for a droplet decreases. For 
T* = 20.7 MeV and ^ = 901 MeV the barrier disappears 
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FIG. 4: Effective potential U{a) as a function of the chiral or- 
der parameter for T = and chemical potential /i = 915 MeV 
(dotted line), fi = 922.7 MeV (sohd line) and fj, = 930 MeV 
(dashed line). 
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FIG. 5: Effective potential U{(j) as a function of the chiral 
order parameter at the critical chemical potential of the first 
order phase transition fi = Hc{T) for temperatures T — 0, 
T = 5 MeV, T = 10 MeV, T = 15 MeV and T = 20 MeV. 



and the line of first order phase transitions ends in a 
critical end point. The two minima merge into one. This 
is completely analoguous to the critical endpoint for the 
water-vapor transition. The computed temperature for 
the endpoint, T^, « 20 MeV agrees well with observation, 
demonstrating the validity of our treatment of the linear 
nucleon-meson model. 

Besides the first order gas-liquid phase transition, the 
effective potential ([s]) also exhibits another first order 
transition at larger values of the chemical potential. Here 
the chiral condensate jumps to much smaller values which 
vanish in the chiral limit ruq = 0. This transition could 
be associated with a transition from nuclear matter to 
quark matter and restoration of chiral symmetry. It is 
obvious that for quark matter the linear nucleon-meson 
model cannot give a valid description and we therefore 
cannot trust Eq. ([8| in this region of the phase diagram. 
It is possible to modify the model in order to incorporate 
an effective change from nucleons to quarks. However, 
the existence and properties of a phase transition depend 
strongly on details of the change of effective degrees of 
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FIG. 6: Chiral order parameter ctq as a function of the chem- 
ical potential for T = (uppermost curve), T = 10 MeV, 
T = 20 MeV, T = 30 MeV, T = 40 MeV, T = 50 MeV, 
T = 60 MeV, r = 70 MeV and T = 80 MeV (lowermost 
curve). 
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FIG. 7: Baryon number density as a function of the chemical 
potential for T = (lowermost curve), T = 10 MeV, T = 
20 MeV, T = 30 MeV, T = 40 MeV, T = 50 MeV, T = 
60 MeV, T = 70 MeV and T = 80 MeV (uppermost curve). 



freedom [TT] such that no reliable information can be 
gained without a better understanding how nucleons are 
replaced by quarks. In practice, the chiral condensate in 
the vicinity of the gas-liquid transition is typically in the 
range a = 65... 93 MeV, while the additional minimum at 
larger chemical potential occurs at much smaller values, 
CT < 10 MeV. Obviously, the Taylor expansion of Uvac for 
a around f^^ in Eq. (|9| is no longer reliable for such small 
values of a. The investigations in ref. ill show that the 
phase transition to quark matter (if it exists) occurs for 
baryon densities higher than the ones that would result 
from the linear nucleon-meson model. Our limitation to 
baryon densities smaller than 1.5 times nuclear density, 
as indicated in Fig. [T] is therefore a conservative estimate 
of the validity of our model. 

For temperatures higher than the one of the critical 
endpoint = 20.7 MeV one finds that the first order 
phase transition gets replaced by a crossover which gets 
rapidly rather smooth. To illustrate this we plot in Fig. [6] 
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the chiral order parameter ctq as a function of the baryon 
chemical potential fi for various values of the tempera- 
ture in the range T = ... 80 MeV. Similarly, Fig. [7] 
shows the baryon density as a function of the chem- 
ical potential for the same temperatures. Together with 
Fig. [2] and Fig. [3] this demonstrates clearly that the ob- 
served chemical freeze-out points in the /i — T-diagram 
are far away from any phase transition or rapid crossover. 
This main result of the present letter contrasts with the 
properties of chemical freeze-out at low baryon density. 

Still, the baryon density (or particle density if mesons 
are included) changes rapidly as a function of tempera- 
ture in the range of interest, cf. fig. [2] The rates of pro- 
cesses with a change of particle numbers, as annihilation 
or production of strange particles, depends strongly on 
the number density of hadrons ^1 ■ These processes stop 
effectively once the hadron density drops below a crit- 
ical value. Since a rather small change of temperature 
corresponds to a substantial change in hadron density 
the chemical freeze-out occurs for a small temperature 
interval. This explains why common values of fi and T 
can describe rather well the abundancies of all hadron 
species. Even if the freeze-out number density for differ- 
ent hadrons differs to some extent, this will not have a 
large effect on the value of the freeze-out temperature. 

We plot in Fig. [T]the line of constant baryon density, 
n = 0.15 nnucioar- The observed freeze-out temperatures 
are well described by this line if the baryon density is high 
enough. We observe that any line with n = crinucioarj 
with c < 1, will end for low T in the line of the gas- 
liquid nuclear phase transition. This simply follows from 
the jump of n from n = to n = n^uclcs^T for T = 0. 
This explains why the observed freeze-out points tend to 
approach the first order nuclear phase transition. As we 
have seen, however, this should not be interpreted as a 
continuation of this first order line by a strong crossover. 

We note as an aside that a freeze-out at constant 
baryon number density can explain also a puzzling fea- 
ture of the chemical freeze-out volume per unit rapidity 
dV/dy. While this quantity grows with collision energy 
for the high energy experiments, the relation is reverse 
for y/s~NN ^ 4.5 Gev [2]. However, assuming constant 
f^Baryons at freczc-out, thc bchavior at AGS energies re- 
fiects just the experimental finding for the energy depen- 
dence of the proton yield dN/ dy at y = [15] . 

In principle, a deviation of the chiral condensate (Tq 
from its vacuum value f-^ and a non-vanishing expecta- 
tion value Wo (implying fXcS = fJ'+g'^a 7^ /^) should lead to 
modified particle yields with respect to a thermal model 
that assumes vacuum masses and /i = /ioff- The chemi- 
cal potential extracted from the thermal fits corresponds 
in first approximation actually to ^ -I- gujQ — h{ao — f-n)- 
However, we find the deviation of this quantity from /i 
to be small in the relevant regime, typically ~ 10 MeV. 

We may also explore the phase diagram for tempera- 
tures higher than the freeze-out temperature. The tem- 
perature dependence of the number density and energy 
density for various values of /i is shown in figs. [8]and[9j 
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FIG. 8: Number density of baryons and pions (solid lines) as 
well as baryons only (dashed lines) as a function of temper- 
ature for the chemical potentials ^ = 550 MeV (lowermost 
curves), = 650 MeV, /i = 750 MeV and = 850 MeV 
(uppermost curves). 

Meson fluctuations 

For our quantitative discussion we have not included 
the fluctuations of the tt, ct, and lo mesons. We also did 
not consider their contribution to the particle density for 
hadrons. For an estimate of the validity of our approx- 
imation we may compute the contribution of the meson 
fluctuations to the effective potential, using a Gaussian 
approximation similar to the baryon fluctuations. For 
that purpose we assume that the momentum dependent 
parts of their inverse propagators are not modified by 
the effects of non-zero chemical potential and tempera- 
ture, and we use a Gaussian or quasi-particle approxima- 
tion. The effect of the bosonic fluctuations can then be 
parametrized in terms of the pressure of a free gas of rel- 
ativistic bosons, Pbg{T, fJ-,m), completely analogous to 
Eq. 0. The masses of the pions and cr-mesons depend 
on a according to 

^'"^ ^a[da^ "''^'^ j ' "^'^"^ = 9a^' ^^^^ 

while TO^ is independent of a in our approximation. 
Thus the w-fluctuations contribute to U only a tem- 
perature dependent constant and we will neglect them. 
The Gaussian approximation becomes invalid in regions 
where iri^{a) or m^(o') become negative. Close to the 
value where m'^{a) vanishes it also matters if one eval- 
uates dU/da for the potential ([s]), as we do it here, or 
uses a self-consistent formulation where U includes the 
meson fluctuations. We stay away from such problem- 
atic regions. 

The pion- and cr-contributions to the effective potential 
A[/„i read 

AC/„ = -3Pbg (T, 0, m,[a)) - Pbg (T, 0, m,(a)) . (16) 

Here Pbg obtains from Eq. ([t]) by setting /i = 0, drop- 
ping the second term and replacing the fermionic by the 
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FIG. 9: Energy density of baryons and pions (solid lines) as 
well as baryons only (dashed lines) as a function of temper- 
ature for the chemical potentials /j, = 550 MeV (lowermost 
curves), fi = 650 MeV, fi = 750 MeV and /i = 850 MeV 
(uppermost curves). 

bosonic mean occupation number (changing +1 to —1 in 
the denominator). We have verified that the meson fluc- 



tuations play only a minor role in the region of interest, 
as shown in Fig. [l]and[2j The relative importance of the 
pion fluctuations can be judged from figs. [8]and[9j 

We conclude that the central result of this note seems 
to be rather robust. The chemical freeze-out at high 
baryon density is not related to any phase transition or 
rapid crossover. It rather follows a simple line of constant 
freeze-out density. We believe that the linear nucleon- 
meson model can give a reliable description for the whole 
low temperature region of the QCD-phase diagram up to 
densities of at least 1.5 times nuclear density. This re- 
quires the inclusion of the meson fluctuations, for exam- 
ple by a genuine functional renormalization group study. 
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